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Abstract 

First-principles density-functional theory calculations show switching magnetiza- 
tion by 90°can be achieved in ultrathin BFO film by applying external electric- field. 
Up-spin carriers appear to the surface with positive field while down-spin ones to 
the negative field surface, arising from the redistribution of Fe-t2# orbital. The half- 
metallic behavior of Fe-3d states in the surface of R phase film makes it a promising 
candidate for AFM/FM bilayer heterostructure possessing electric- field tunable FM 
magnetization reversal and opens a new way towards designing spintronic multifer- 
roics. The interface exchange-bias effect in this BFO/FM bilayer is mainly driven 
by the Fe-t2# orbital reconstruction, as well as spin transferring and rearrangement. 

Key words: Magnetoelectric coupling; Ultrathin BFO film ; Spin transfer ; 
Antiferromagnetic structure 
PACS: 73.61.-r,75.70.-i,71.15.Mb 



* Corresponding author. 

Email address: hjfeng@nwu.edu.cn, fenghongjian@gmail.com (Hong- Jian 

Feng) . 



Preprint submitted to Elsevier 



1 Introduction 



BiFe03(BF0) attracted much attention due to the coexistence of antiferro- 
magnetic(AFM) and ferroelectric(FE) behavior, as well as the enhanced weak 
ferromagnetic (FM) magnetization and FE polarization in films on SrTiOs sub- 
strates prepared by pulsed laser deposition(PLD)[I]. The enhanced FM and 
FE properties in BFO films are attributed to the suppression of cycloidal 
magnetic structure with a long wavelength period of ~ 620 Aj2] caused by 
the Dzyaloshinskii-Moriya interaction(DMI) in bulk pJlH] . DMI is able to pro- 
duce the canting of neighboring Fe ions with G-type AFM structure observed 
under a Neel temperature of ~ 643 K[5], along with the spiral magnetiza- 
tion along the crystal. In terms of magnetization enhancement, DMI gives 
rise to the tilting of the neighboring G-type AFM Fe ions leading to the 
weak FM magnetization, and the side-effect is the spiral spin structure which 
tends to reduce the macroscopic magnetization. Fortunately, doping [6. 7, 8j and 
thin film fabrication pQ can be used to impede the spiral magnetization owing 
to the suppression of magnetic cycloid by doping ions and size-effect in epi- 
taxial films. However, from a practical point of view several problems have 
to be solved: (i) The weak magnetoelectric coupling is quadratic due to the 
spin cycloid in the single phase BFO pfTUlfTT] . (ii)The indirect magnetoelectric 
coupling is mediated by the direct coupling between antiferrodistortive(AFD) 
distortions being the out of phase rotation of neighboring oxygen octahedra 
along [1 1 1] direction and AFM vectors fT2lll3lll4j . (hi) Electric-field controlled 
FM magnetization is not able to be implemented because of the FE-AFM 
coupling in single phase BFO. In order to drive FM magnetization rever- 
sal by external electric-field in magnetic memory storage with lower energy 
consumption, a heterostructure composed by the BFO AFM layer and FM 
layer, e.g. BFO/La .7Sr .3MnO(LSMO) heterostructure [T5l[T6] .is proposed to 



3 



accomplish it instead of using the single phase BFO films where the quadratic 
magnetoelectric coupling is FE-AFM interaction |T7TfT8] . Electric-field tunable 
FM magnetization in BFO/FM heterostructure is mainly achieved through 
exchange-bias(EB) effect in the heterointerface [T9ll20ll21j .and the EB coupling 
effect is more than orbital reconstruction, the AFM magnetization rotation, 
spin transfer, and electronic structure in the interface are all playing a role 
to it, so exploring the above-mentioned properties in BFO film would be a 
significant step towards understanding the origin of it. Motivated by the skin 
layer recently observed in the BFO single crystals [22]. we proposed an ultrathin 
BFO film with considering the rhombohedral R3c (R) phase and the tetragonal 
P4mm(P) phase recently observed in the morphotropic phase boundary(MPB) 
in BFO film[23j. We carried out density functional theory(DFT) based calcu- 
lations to investigate the AFM magnetization rotation and the microscopic 
mechanism involved under external elect ric-field. 

2 Computational details 

We use the local spin density approximation(LSDA) to DFT scheme with 
a uniform energy cutoff of 500 eV as in our previous work [T2ll 1311 141123] . We 
considered Bi 5s, 5p, and 6s electrons, Fe 3s, 3p, and 3d electrons, and O 2s and 
2p electrons as valence states and use llxllxl Monkhorst-Pack sampling of 
the Brillouin zone for all calculations. We used the slab model to construct (0 
1) surface in hexagonal frame of reference for R and P phase film having same 
vacuum region as the atomic layers. The relaxation was carried out with the 
top and bottom three atomic layers being moved as the forces on the ions were 
less than 0.01 meV/A. G-type AFM spin configuration was used to construct 
the initial AFM structure, and the angel between the AFM direction and the [1 
1 1] FE polarization direction in rhombohedral representation was changing 
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from 0°to 180°to calculate the magnetic anisotropic energy(MAE) and the 
favorable spin configuration with considering the spin-orbit coupling(SOC) 
and noncollinear magnetism. In addition, considering the FM ordering with 
the same condition as mentioned in AFM case, we set the same variation 
of angle between FM spin and [1 1 1] direction to calculate the MAE in 
comparison with AFM structure. We did not introduce the on-site Coulomb 
interaction by adding a Hubbard-like term to the effective potential in that 
DFT+U scheme is susceptible to the DMI[12], and further the MAE. A saw- 
tooth like external electric potential f25ll26] has been used as 

V ext (r) = 47rm(r/r m - 1/2), < r < r m (1) 

where m is the surface dipole density of the slab, r m is the periodic length 
along the direction perpendicular to the slab. We used electric-field of 2 V/A 
for R phase film and 1 V/A for T phase film due to the restoring forces and 
overshoot in the iteration process as electric-field is beyond that value. 

3 Results and discussion 

The initial structure for R and P phase films are taken from Ref. [1], then 
ions in the vicinity of the surface are all fully relaxed under field. The crys- 
tal structure of the R phase ultrathin film before and after applying field is 
shown in Fig.l. It is apparent that the Bi ions are displacing significantly 
along the field comparing with the case in Fe, implying the Fe ions are con- 
fined in oxygen octahedra and slightly affected by the field. Interestingly, with 
applying an elect ric-field, Fe ions in P phase films are actively participating 
in the tetragonal distortions owing to the epitaxial strain in-plane. AFM or- 
dering is energetically favored for R phase films before and after exerting a 
field ,as well as for P phase films under elect ric-field, but FM ordering is more 
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favorable in P phase films without external bias due to the unreleased strain 
in-plane. Considering about the anisotropic superexchange interaction [273 an d 
single-ion anisotropy[28], the MAE per formula unit before and after applying 
external bias were reported for R and P phase in Fig. 2 and Fig. 3, respectively. 
Sinusoidal waveform of MAE is established for AFM spins except the P phase 
films with FM ordering where anomalies exist around the peak showing that 
the FM spins are not stable comparing with the AFM ones, and the gaus- 
sian fitting curves are demonstrated in Fig. 4 for this case. Without applying 
bias, the R phase films exhibit an AFM spins parallel to the [1 1 1] polar- 
ized direction in rhombohedral representation, then spins rotate to direction 
perpendicular to the [1 1 1] direction under field. Moreover, the FM spins 
parallel to the [1 1 1] direction is also rotating to the direction perpendicular 
to the [1 1 1] in the P phase films with applying electric- field. Eventually, an 
electric- field driven (111) easy magnetization plane is established in these 
two phase films, implying that a rotation of magnetization by 90°is able to be 
induced by the bias along [1 1 1] direction, and this phenomenon is expected 
to be observed in the skin layer of BFO bulk. From Fig. 4 the gaussian fitting 
curves indicate that a local minimum MAE is calculated at ~ 60°relative to 
the [1 1 1] direction, and this FM ordering is a metastable sate and prone to 
be overcome by the electric- field. The electric-field tuned rotation of magneti- 
zation and spin rearrangement in the ultrathin BFO film is depicted in Fig. 5. 
A heterostructure composed by this ultrathin BFO film and thin FM film 
is expected to possess an electric-field controlled FM magnetization reversal 
through the interface EB coupling. 

In order to demonstrate further the magnetization and spin transfer behavior 
across the films under external bias, we reported the magnetization averaged 
in-plane across the R phase and P phase films in Fig. 6 and Fig. 7, respectively. 
One can see that the up-spin and down-spin carriers are equally spaced across 
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the films, and they are equal in magnitude and opposite in direction which is 
caused by the AFM spin configuration. The maximum spin density is about 
0.2 11 b in R phase films which is less than 0.35 /x^ in P phase. The relative large 
value of magnetization in P phase is attributed to the epitaxial strain in-plane, 
and the tetragonal displacement of Fe ions out-of-plane which tends to couple 
with the oxygen octahedra rotation to produce a much higher spin density. 
We define the direction of a surface as the normal pointing out from it and the 
sign of the field as the cosine of the angle between the direction of the field and 
surface, then there must be a positive field and negative field in two sides of 
a film for a given field. From Fig. 6 and Fig. 7, it is worth mentioning that the 
up-spin carriers are always moving along the electric-field while the down-spin 
carriers are in the opposite direction. Eventually, the up-spin carriers appear 
in the surface where positive field exists, leaving the down-spin ones in the 
other side of the film with negative field. When the BFO skin layer is combined 
with the FM layer to achieve the bilayer heterostructure, the up-spin carriers 
are accumulating in the surface with positive electric-field as the down-spin 
ones aggregating in the opposite side. The spin carriers are coupling to the FM 
layer through EB involving the orbital reconstruction and the spin transferring 
behavior, hereby the spin carriers would be reversed by the alternating electric- 
field leading to the electric-field driven FM magnetization reversal. 

To shed light on the surface spin transfer and intrinsic electronic property, the 
density of states(DOS) of triply degenerate Fe-t 2 ^ orbital in R and P phase 
film surface with negative field are illustrated in Fig. 8 and Fig. 9, respectively. 
In BFO the Fe-3d are decomposed into triply degenerate t 2g and doubly degen- 
erate e g states due to the crystal field splitting, and the latter is hardly to be 
find in the surface showing that the spin carriers are mostly arising from the 
former. It is clear that the occupied Fe-t 2 ^ states are changing from majority 
spins to minority ones, indicating that the down-spin carriers are accumulat- 
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ing in the surface with negative field , and the corresponding phenomenon is 
able to be confirmed in the surface with positive field. The accumulation and 
rearrangement of spin carriers is determined by the redistribution of Fe-t2 9 
orbital occupation with exerting external bias. Moreover, there exists a small 
amount of minority spins found both in R and P phase films without external 
bias, and the majority spins are still dominating in the film. On the contrary, 
all the occupying states are minority spins under electric-field, this partly re- 
flects the accumulating behavior in the surface which is consistent with the 
spin density calculation shown in Fig. 6 and Fig. 7. It is worth pointing that 
there are finite DOS of up-spin and down-spin Fe-t 2g states passing through 
the fermi level in R phase film before and after turning on the bias, clearly 
showing the half-metallic property in Fe-3rf states. The half-metallic behav- 
ior achieved in the ultrathin BFO film enables it a potential candidate for 
spintronic applications. Without field the occupying states are coming from 
up-spins in R phase film while the up-spin states passing through fermi energy 
is able to provide a transferring way across the surface, and the transferring 
path is provided by the down-spin ones under field. Therefore we suggest the 
half-metallic property in R phase film constitutes the interface EB coupling 
mechanism, and the R phase ultrathin BFO film is proposed to be an ideal 
candidate for making electric-field tunable FM magnetization reversal devices 
and spintronics. 

4 Conclusion 

DFT based first-principles calculations have been performed to study the mag- 
netization rotation and spin density in ultrathin BFO film. A (1 1 1) easy mag- 
netization plane is achieved for R and P phase films with applying external 
bias, hereby an AFM magnetization rotation driven by electric-field is accom- 
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plished. The up-spin carriers are along the field direction and accumulating in 
the surface with positive field while the down-spin ones are opposing the field 
and aggregating on the other side of the film. The half-metallic Fe-3c? states 
make R phase BFO ultrathin film a prosperous candidate for spintronics. 
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Figure Captions: 

Fig.l The structure of the R phase BFO ultrathin film before and after 
exerting field. 

Fig. 2 Dependence of total energy on the angle between spins and [1 1 1] 
polarized direction in R phase film with and without applying external bias, 
respectively. 

Fig. 3 Dependence of total energy on the angle between spins and [1 1 1] 
polarized direction in P phase film with and without applying external bias, 
respectively 

Fig. 4 Gaussian fitting of MAE in P phase film without applying electric-field. 

Fig. 5 Schematic diagram for the magnetization rotation under electric-field, 
the dotted line denote the [1 1 1] direction, and the field is applying along 
the [1 1 1] direction. 

Fig. 6 Averaged spin density in-plane across the R phase ultrathin BFO film 
under external bias. 

Fig. 7 Averaged spin density in-plane across the P phase ultrathin BFO film 
under external bias. 

Fig. 8 t 2g DOS for R phase film before and after applying electric-field. 
Fig. 9 t2g DOS for P phase film before and after applying electric-field. 
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Fig. 1. The structure of the R phase BFO ultrathin film before and after exerting 
field. 
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Fig. 2. Dependence of total energy on the angle between spins and [1 1 1] polarized 
direction in R phase film with and without applying external bias, respectively. 
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Fig. 3. Dependence of total energy on the angle between spins and [1 1 1] polarized 
direction in P phase film with and without applying external bias, respectively. 
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4. Gaussian fitting of MAE in P phase film without applying electric- field. 
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Fig. 5. Schematic diagram for the magnetization rotation under electric- field, 
dotted line denote the [1 1 1] direction, and the field is applying along the [1 
direction. 
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Fig. 6. Averaged spin density in-plane across the R phase ultrathin BFO film under 
external bias. 
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Fig. 7. Averaged spin density in-plane across the P phase ultrathin BFO film under 
external bias. 



17 




8. t2 9 DOS for R phase film before and after applying electric- field. 




9. t2 9 DOS for P phase film before and after applying electric- field. 
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